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Abstract 



J>v The properties of the most probable ground state candidate for the XXZ spin 

u , chain with the anisotropy parameter equal to —1/2 and an odd number of sites 

is considered. Some linear combinations of the components of the considered 
CN ■ state, divided by the maximal component, coincide with the elementary symmet- 

ric polynomials in the corresponding Bethe roots. It is proved that those polynomi- 
| als are equal to the numbers providing the refined enumeration of the alternating- 

' sign matrices of order M + 1 divided by the total number of the alternating-sign 

matrices of order M, for the chain of length 2M + 1. 

o : 

in ■ 

§ '■ 1 Introduction 

° ■ 

The XXZ spin chain is a system consisting of N spin one-half particles arranged in 
a linear structure. Here each particle interacts with its nearest neighbors only. The 

^ ■ state space V of the XXZ model is the tensor product of N copies of the state space of 
a single spin one-half particle, which is isomorphic to the complex two dimensional 
coordinate space C 2 . Denote the vectors of the standard basis of the vector space C 2 

j> ' as |+) and |— ), so that 

f '->-(?: 

The standard basis of the state space of the XXZ model is formed by the vectors 

1^1,^2/- ..,Pn) = ® ® •••® 

where the variables ■ ■ ■ , }In take the values + and — . Sometimes it is convenient 

to assume that such variables take the values +1 and —1. The used interpretation is 
always clear from the context. 

Denote by cr x , cr^ and cr z the usual Pauli matrices and the corresponding operators 
in C 2 . For any n = 1,2, ... ,N define 

cr l n = I g • • • g> I ® cr 1 ® I . . . g L i = x,y,z. 
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Here we denote by I the identity operator in C 2 . Now we are able to write the Hamil- 
tonian of the XXZ spin chain. It has the form 

-i N 

Hxxz(A) = --£ [crX+i + °%o* +1 . (!) 

;=i 

where Z\ is the anisotropy parameter. In the present paper we use the periodic bound- 
ary conditions. 

As for any quantum mechanical system, the main task related to XXZ spin chain 
is to find the eigenvalues and eigenvectors of the Hamiltonian. Before discussing the 
methods used to find them, let us make two general remarks. Note that the operator 

N 
n=l 

commutes with Hxxz(^)- Therefore, one can look for the eigenvectors of the Hamilto- 
nian Hxxz(^)/ restricting oneself to a subspace spanned by the eigenvectors of E with 
a fixed eigenvalue. It is clear that 

L\ji lt ...,ji N ) = (^Hn \ If1'--->/*n) = (N+ -N-)\iii,...,ji N ), 

where N + and N are the numbers of pluses and minuses,respectively, in the set 
formed by . . . , Hence we can look for the eigenvectors of the Hamiltonian 
Hxxz(^) in a subspace spanned by the basis vectors with a fixed number of pluses, or 
minuses. Considering such subspaces, we will denote by K the number of minuses. 

Note here that the minimal eigenvalue in a subspace with a fixed K is non- 
degenerate [1 . Theorem 4]. 

It is also not difficult to get convinced that the Hamiltonian Hxxz(A) commutes 
with the operator 

N 

p = EK- 

n=l 

One easily obtains that 

P • • • , Pn) = \-V-ii- ■ ■/— ^n)- 

Therefore, each eigenvector \Y) in the subspace with K minuses has the companion 
P \ Y) in the subspace spanned by the basis vectors with N — K minuses. The both 
eigenvectors certainly correspond to the same eigenvalue. In the case of an even N = 
2M and K = M two eigenvectors connected by the operator P belong to the same 
subspace and can coincide. 

It is clear that the subspaces with K = and K = N are one dimensional. Therefore, 
the corresponding basis vectors must be eigenvectors of Hxxz(^)- Thus, we know at 
least two eigenvectors of Hxxz(^)- These eigenvectors are ground states for the XXZ 
spin chain only for the case A > 1. They do not depend on A, and the corresponding 
correlation functions are not quite interesting from the physical point of view. 

For A < 1 and an even N = 2M the ground state of the XXZ spin chain is unique 
and belongs to the subspace with K = M, see, for example, the monographs (2H3|. For 
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A < 1 and an odd N = 2M + 1 we have two ground states connected by the operator 
P, one in the sector with K = M, and the other one in the subspace with K = M + 1. 

In general the eigenvectors and eigenvalues of the Hamiltonian Hxxz(^) can be 
found with the help of the famous Bethe ansatz |4J. Unfortunately, this method pro- 
vides us with the expressions of the eigenvectors and eigenvalues via the roots of the 
Bethe equations, which are not known usually 

Another method, partially equivalent to the Bethe ansatz, is based on the relation 
of the Hamiltonian Hxxz(^) an d the transfer-matrix of the statistical six- vertex model. 
The eigenvalues of the transfer matrix satisfy the famous scalar Baxter T-Q equation, 
see Chapters 9, 10] and |5|. As was proved in paper |6|, in the case A = —1/2 and 
an odd N = 2M + 1 this equation has a solution corresponding to an eigenvalue of 
the Hamiltonian Hxxz(^) equal to — 3N/4. A remarkable generalization to the case of 
the eight- vertex model when it describes the off-critical deformation of the A = —1/2 
six- vertex model is under active study now, see J7H8H2I and references therein. 

It is very likely that the solution of the scalar Baxter T-Q equation found in the 
paper |6[ corresponds to the ground state of the system. Using MATHEMATICA, we 
found explicitly the eigenvectors with the eigenvalue — 3N/4 for N < 17. Many re- 
markable relations between the components of the obtained eigenvectors were noticed 
and formulated in our paper |10| as conjectures. It turned out, for example, that if we 
normalize the eigenvector so that the minimal component is 1, then the largest com- 
ponent is equal to the number A(M) of the alternating-sign matrices of order M. We 
refer the reader to the book ifTTTl and paper [12J for information on the alternating-sign 
matrices. Some of our conjectures are already proved fl3llT4i ri5l, and some of them 
are generalized to the case of different boundary conditions I161 IT7I . 

In the present paper we continue our study of the components of the ground state 
candidate for the XXZ spin chain with A = —1/2 and an odd N. In Section 2 we ob- 
serve a possible connection of some sums of components of the ground state candidate 
with the numbers describing the refined enumeration of the alternating-sign matrices. 
In Section 3 we relate those sums with the elementary symmetric polynomials in the 
Bethe roots, and in Section 4 we prove that these symmetric functions are related to 
the refined enumeration of the alternating-sign matrices. 

2 Empirical results 

In this section we consider the explicit form of the components of the ground state 
candidate of the Hamiltonian Hxxz(^) f° r ^ = —1/2 and an odd N = 2M + 1. 
Note that the shift operator S, defined by the equality 

S |^l,^2/ ••• / Pn-v Hn) — \V-2i ■ ■ ■ / ^N-l/ V-n, Mi)/ 
and the reflection operator R, defined by the equality 

R . . ., pn-i, Pn) = \y-Ni V-tss-1/ ■ ■ ■ ' Y-ii V\)> 

commute with the Hamiltonian Hxxz(^)/ an d leave the subspaces spanned by basis 
vectors with a fixed number of minuses invariant. As we already noted, for the case 
A = — 1/2 an odd N = 2M + 1 there are two ground states, one in the subspace with 
K = M, and another one in the subspace with K = M + 1. They are non-degenerate 
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in the corresponding subspaces. Hence, each ground state is an eigenvector of the 
operators S and R. Actually they are invariant with respect to the action of S and R 
that results in equality of the corresponding components. The ground state candidates 
under consideration are invariant with respect to the action of the operators S and R 
by construction. 

Consider the ground state candidate \Y) belonging to the subspace with K = M. 
It is convenient to redenote the corresponding basis vectors as \n\, rii,..., Hm), where 
the integers n\,ni, . . .,Hm satisfy the condition 1 < n\ < ni < . . . < < N and 
describe the positions of the minuses in the initial notation for the basis vector under 
consideration. The components \Y) are also marked by the integers n\, rii,..., n^. 

Consider the case N = 11. In this case the independent components of the ground 
state candidate are [10] 



ifl,2,3,4,5 


= 1, 


if 1,3,4,5,6 


= 5, 


1^1,2,4,5,6 


= 10, 


ip-1,4,5,6,7 


= 11, 


if 1,3,5,6,7 


= 34, 


if 1,3,4,6,7 


= 41, 


ipl,3,4,5,7 


= 23, 


vp-1,2,5,6,7 


= 30, 


if 1,2,4,6,7 


= 60, 


if 1,5,6,7,8 


= 14, 


ipl,4,6,7,8 


= 52, 


url,4,5,7,8 


= 73, 


if 1,4,5,6,8 


= 46, 


if 1,3,6,7,8 


= 75, 


ipl,3,5,7,8 


= 169, 


if 1,3,5,6,8 


= 128, 


if 1,3,4,7,8 


= 101, 


if 1,2,6,7,8 


= 42, 


ipl,2,5,7,8 


= 114, 


url,4,7,8,9 


= 81, 


if 1,4,6,8,9 


= 203, 


if 1,4,6,7,9 


= 174, 


1^1,4,5,8,9 


= 141, 


if 1,3,6,8,9 


= 226, 






if 1,3,5,8,9 


= 260, 


uyl,3,5,7,9 


= 429. 







Here we normalize the vector in such a way that the minimal component Y 1,2,3,4: ' 5 
is equal to 1. The components not given above can be obtained using the shift and 
reflection invariance. 

In accordance with the conjecture formulated in the paper |10| the maximal com- 
ponent Y 1 ' 3 ' 5 ' 7 ' 9 i s equal to the number A (5) of the alternating-sign matrices of order 
5. But this is only the tip of an iceberg. 

Let us consider (?) = 5 components which can be obtained from the maximal 
component Y 1,3,5,7,9 by incrementing one of its indices by 1. We have 

if 2,3,5,7,9 = 169/ if 1,4,5,7,9 = ^ if 1,3,6,7,9 = ^ 
if 1,3,5,8,9 260 if 1,3,5,7,10 429 

The sum of the five obtained numbers is equal to 1287. 

It is the time to look at Figure |l] borrowed from the paper by Bressoud and 
Propp |12J. Let us explain its meaning. 

Recall that an alternating-sign matrix is a matrix with entries 1, 0, and —1 such that 
the 1 and —1 entries alternate in each column and each row and such that the first and 
last nonzero entries in each row and column are 1. One can easily get convinced that 
the first column of an alternating-sign matrix contains only one entry 1 and all other 
entries are zero. Denote by A(M,r) the number of alternating-sign matrices of order 
M for which the unique entry 1 is at the rth position in the first column. These are the 
numbers that are given in the rows of Figure |T] One says that the numbers A(M,r) 
give the refined enumeration of the alternating-sign matrices. 
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1 

1 1 

2 3 2 

7 14 14 7 

42 105 135 105 42 

429 1287 2002 2002 1287 429 

7436 26026 47320 56784 47320 26026 7436 

Figure 1: The numbers A(M, r) 

Pay our attention to the sixth row of Figure [U We meet there the numbers 429 
and 1287 already familiar to us. Let us go further and consider the components of the 
ground state which can be obtained from the component W 1 > 3 > 5 > 7 > 9 by incrementing two 
its distinct indices by 1. The possible (2) = 10 variants give 



^2,4,5,7,9 


= 128, 


^2,3,6,7,9 


101, 


^2,3,5,8,9 _ 


= 114, 


^2,3,5,7,10 


= 203, 


^1,4,6,7,9 


= 174, 


u/4,4,5,8,9 __ 


141, 


^1,4,5,7,10 


= 226, 


^1,3,6,8,9 


= 226, 






^1,3,6,7,10 _ 


260, 


^1,3,5,7,10 


= 429. 







The sum of the obtained components is equal to 2002. We again obtain a number from 
the fifth row of Figure HI 

Taking into account the shift and reflection invariances, one can see that (|) = 10 
components which can be obtained by incrementing three of the five indices are nu- 
merically coincide with the corresponding components obtained by incrementing two 
indices, and (f) = 5 components which can be obtained by incrementing four of the 
five indices are numerically coincide with the components obtained by incrementing 
one index. Furthermore, decrementing indices we obtain the same set of numbers as 
incrementing them. 



3 Symmetric polynomials in Bethe roots 

Let Sk be the symmetric group of degree K. According to the Bethe ansatz |4J the 
eigenvectors of the Hamiltonian Hxxz(^) f° r general N and K have the form 

W= E YHl ' n2 nK \ni,n 2 ,...,n K ) 

l<ni<n 2 <---<n K <N 

L E A B z% 1 f% 2y ..z n s fc\n 1 ,n2,...,n K ). (2) 

l<tii<n 2 <--<n K <N seSjc 

Here z\, Zi, . . ., z& are complex numbers, satisfying the Bethe ansatz equations 

Z N = { _ 1)K -i Yl [f( Zl>Zk )/f( Zk>Zl) ] , k = l,...,K, (3) 
1=1 

with 

/(zi,z 2 ) = 1 - 2Az 2 + ziz 2 . 
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We call the numbers Z\,zi, . . .,Zk, satisfying ©, the Bethe roots. The coefficients A s 
have the form 

A s = sgn(s) Yl f( z s(h)' z s(h))' 

l<k 1 <k 2 <K 

where sgn(s) is the sign of the permutation s. The eigenvalue of Hxxz(^)/ correspond- 
ing to the eigenvector (01, is 

E = -±AN+f^(2A-z k -z k 1 ). (4) 
z k=l 

It is convenient for our purposes to introduce the notation 

T>ni,n 2 ,...,n K _ A n x n 2 „n K 

° s ~ /lsZ s(i) z s(2) ■ ■ -^(rr 

Using this notation, we can write 

seS K 

Return again to the case of an odd N. Consider an eigenvector \W) belonging to 
the subspace with K = M, not necessary a ground state. Denote by F the component 

Y 1 ' 3 2M -\ so that 

f = £ bp 2M -\ 

seS M 

Incrementing an index of Y 1 ' 3 '-' 2M ~ 1 by 1 we obtain M components 

r _ xp-1,3,— ,2m-3 7 2m,2m+l,...,2M-l 

It is clear that 

p _ V -1 7 Dl / 3 / ...,2m-3,2m-l,2m+l,...,2M-l 
t m — ^s(m) D s i 

seS M 

and after summation over m we have 

M M 

E c V V -r Dl,3,—/2Jw-3 / 2»i-l,2m+l,...,2M-l 
= 2^ 2^ Z s(m)^s 
m=l m=lseSM 

/ M \ 

= £ £ Z gl,3 / ...,2jn-3,2OT-l,2m+l / ...,2M-l 



M \ 

y z j y gl,3,---/2m-3,2m-l / 2m+l,...,2M-l 
v m=l / seS M 



where we use the evident equality 



M M 

z s(m) = Zm ' 
m=l m=l 

Thus, for the first symmetric polynomial |18| in the Bethe roots we have the expression 

M i M 

e\(z\,. . . ,zm) = YL Zm = pJl Fm - 

m=l k=\ 



Denote now by 7 mx ,...,m r the component which is obtained from *plA— /2M-1 j^y 
incrementing the indices at the positions m\,...,m r by 1. Here we assume that 
1 < m\ < . . . < m r < M. In a similar way as above, using the equality 

z s(m 1 ) ■ ■ ■ z s(m 1 ) = ' ' ' Zm r' 

l<m\<...<m r <M \<m\<...<m r <M 

we see that the rth elementary symmetric polynomial in the Bethe roots is related to 
the components of the eigenvector \W) as 

&r{z\r • • • /%m) = 5^ Z mi Z m2 . . . Z mr = — ^m\,mi,...,m r - 

l<mi<m2<...<m r <M l<)« 1 <?«2<...</K r <M 

In the next section for the case of A = —1/2 and an odd N we show that for the 
considered state with the eigenvalue — 3N/4 the elementary symmetric polynomials 
in the Bethe roots coincide with the numbers giving the refined enumeration of the 
alternating-sign matrices. 



4 Relation to refined enumeration of alternating-sign 
matrices 

It is well known that the XXZ spin chain is closely related to the six-vertex model, 
see, for example, the monograph [2 J. To describe this relation, let us parametrize the 
Boltzmann weights of the six-vertex model in the following way 

a(u) = (t(tu), b(u) = <t{t~ 1 u), c(u) = ct(t 2 ), 

where t is a fixed parameter, u is a free parameter, and we use the convenient notation 

<j{x) = x — x^ 1 

introduced by Kuperberg. One can show that the Hamiltonian Hxxz(^) is connected 
with the commuting family T(u) of transfer matrices of the six- vertex model by the 
equality 



T~\u) 



dT(u) 
du 



T(T 2 — T 



-2> 



Hxxz(^(t))-^(t 2 + t- 2 ) 



where 



A(t) = (t 2 + t- 2 )/2. 



(5) 



(6) 



Therefore, one can diagonalize the Hamiltonian Hxxz(^) with A given by ® and the 
family T(u) simultaneously 

It appears that the vector \Y), defined by relation (J2J), is an eigenvector of T(u) with 
the eigenvalue 



A(u) = a N (u)U(u,zi) . . . U(u,zk) + b N (u)V{u,z{) . . . V(u,zk), 



N, 



(7) 



where 



U(u,z) 



a(u)b(u) + {c 2 (u) — b 2 (u))z 
a 2 (u) — a(u)b(u)z 



V{u,z) 



a 2 {u) — c 2 {u) — a(u)b(u)z 
a(u)b(u) — b(u) 2 z 



see, for example, the paper [19J and the monograph |2J. It is convenient to introduce 
the complex numbers u kl k = 1, . . . , K, related to the Bethe roots z k by the equality 

a(u k ) a{ru k ) T 2 u 2 k - 1 
k b{u k ) cr(T-l Mfc ) u 2 -t 2 ' W 

It is clear that the numbers are defined up to a sign. One easily obtains 
U{u,z k ) = ; — ^ — , V{u,z k ) = — = s— . 



Denoting 1 



cr(uu k ) " a(uu k 



K 

Z(u)=H[u k (r{uu k 1 )}, (9) 
fc=i 



we rewrite equality as 

A(u)^(u) = ct n (to)^(t- 2 i/) + ct n (t- 1 m)^(t 2 m). (10) 

Treating ^(m) as eigenvalues of a commuting family of the operators Q(u), we come to 
the equality 

T(u)Q(u) = (t n (tu)Q(t- 2 u) + (t n {t- 1 u)Q{t 2 u), 



which is the famous Baxter T-Q equation. The equation (jTOjl is called the scalar Baxter 
T-Q equation. One can solve equation dTol instead of solving the Bethe ansatz equa- 
tions ©. The corresponding Bethe roots z k can be recovered from the roots u k of the 
Laurent polynomial using relation ©. 

Assume now that T = exp(izr/3) that gives A = —1/2. It appears that for an 
odd N = 2M + 1 the Bethe equations in the subspace with K = M has a solution 
corresponding to 

r iN 
A(u) = \a(u) +b(u)] N = (t + t- 1 )(u-u- 1 ) =cr N (u). 



Here we use the identity 



T + T- 1 = 1. 



To find the solution in question, one observes that defining 

<p(u) = cr N (u)£(u) 

and using the equality 

ct(t 3 x) = —a(x), 

one can rewrite equation dTol as 

<p(u) + cp(T 2 u) + (p(rSi) = 0. (11) 

1 Actually the function £(u) is defined up to a factor. We define it in such a way that the coefficient at 
the maximal degree of u is equal to 1 . 
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Let us suppose that the solution under consideration corresponds to a ground state 
of the system in the subspace with K = M. As was proved in the paper |1J, in this case 
one can number the Bethe roots in such a way that 

One can chose the corresponding numbers u m so that 

-l _ 

In this case the function £(«) satisfies the equality 
that gives 

(p{u~ l ) = —<p{u). (13) 

It is also clear that <p(u) should be a centered Laurent polynomial in u of degree 6M + 2 
and that 

<p(-u) = (-l) M+1 cp(u). (14) 

Finally note that the Laurent polynomial <p(u) should be divisible by cr 2M+1 (u) by 
construction. 

It was shown in the paper [6 J that up to a constant factor there is a unique solution 
of equation (JTTJl which is a centered Laurent polynomial in u of degree 6M + 2, satisfies 
relations dl3l and ((14)1 , and is divisible by (7 2M+1 («). This polynomial has the form 

We chose the normalization in accordance with the definition (|9jl of £(«). 

Let us now relate the Laurent polynomial <p(w) to elementary symmetric polyno- 
mials in the Bethe roots. Define the variable 

a(u) <t{tu) t 2 u 2 — 1 



b(u) ct(t _1 w) u 2 — T 2 
then, using relation one obtains 

Cr(T _1 w)<7(TM^ ) 

Consider the function 

M M 



(16) 



x( z ) = n( z_zm ) = E(-i) r e r (zi,. . .,zm)z 

m=l r=0 



M-r 



where we assume that eo(z\, . . . ,Zm) = 1- Using ((16} , we come to the equality 

x(z)-f\( Z -z ) ^> 2M+ VM") a7) 

9 



It follows from the explicit representation ((15)) for the function <p(u) that it satisfies 
the second order differential equation 

u ± (u**p.) - 6M l 4^u^P- + (3M + 1) (3M - l)cp(u) = 0. 
du V du J u — 1 du v ;v ;rv ; 

Now one can show that the function x( z )> related to <p(u) by equality (fT7)l , satisfies the 
differential equation 

z(z + l)^r^ + 2(z - M)^^ - M(M + l)#(z) = 0. 

Using this equation, we come to the recursive relation for the elementary symmetric 
polynomials in the Bethe roots, 

(2M-r + l)re r (z 1/ ...,z M ) = (M-r + l)(M + r)e r „ 1 (z 1/ ...,z M ). 

Comparing it with the recursive relation for the numbers giving the refined enumera- 
tion of the alternating-sign matrices 

(2M-r-l)rA(M,r + l) = (M - r)(M + r - 1) A(M,r), 

see, for example, |20|, and using the initial condition e§{z\, . . ■ ,Zm) = 1/ we see that 

e r ( Zl ,...,ZM) = ™yA(M + l,r + l). 

Hence, if our conjecture, stating that if we normalize the ground state candidate so 
that the minimal component is equal to 1, then the largest component is A(M) IITOl , 
we come to the relation 

E F mi ,m 2 m r = A(M + l,r+l). 

l<mi<m2<:.<tn r <M 

For the simplest symmetric polynomial in the Bethe roots, e\{z\, . . .,Zm), we have 

M 

e 1 (z lr ... / z M ) = £ Zm = (M + l)/2. 

m=\ 

Taking into account that in the case under consideration 

M M 

Ylt Z m = Ylt Zm ' 

m=l m=l 

we obtain from equality l|4jl that the eigenvalue of Hxxz(A), corresponding to the dis- 
cussed ground state candidate, is equal to — 3N/4. Certainly, this result can be ob- 
tained in a more direct way, starting from relation 
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5 Discussion 



In the present paper we considered the properties of the most probable ground state 
candidate for the XXZ spin chain with the anisotropy parameter A equal to —1/2 and 
an odd number of sites N. We related some combinations of the components of the 
considered state to the elementary symmetric polynomials in the corresponding Bethe 
roots. Then we proved that those polynomials coincide with the numbers giving the 
refined enumeration of the alternating-sign matrices. 

The components Y ni, - ,nK , defined via equality (|2jl are antisymmetric polynomials 
in the Bethe roots z\, . . . , z%. Representing them as 

y mi nK =( n (^-zgW nK , 

\l<k 1 <k 2 <K J 

we introduce the quantities <p n i'—> n z which are symmetric polynomials in the Bethe 
roots. Therefore, the components of an eigenvector of the Hamiltonian Hxxz/ given by 
the Bethe ansatz (0, can be expressed, in principle, via elementary symmetric poly- 
nomials in the Bethe roots. Thus, the components of our ground state candidate can 
be expressed via the numbers giving the refined enumeration of the alternating-sign 
matrices. 

The authors of the paper [21\, starting from relation (fl5)> , found, in particular, a 
representation of the elementary symmetric polynomials in the Bethe roots via some 
sums of products of binomial coefficients. 

Special solutions to the Baxter T-Q equation can be found for different boundary 
conditions [22 , 23J, and our results can be generalized to those cases. 
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